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Abstract. - Depletion interactions between colloidal particles dispersed in a fluid medium 
are effective interactions induced by the presence of other types of colloid. They are not 
instantaneous but build up in time. We show by means Brownian dynamics simulations that the 
fluctuations in the depletion force between two guest particles in a host dispersion of differently 
sized colloids do not decay exponentially with time, but show a power-law dependence. A simple 
scaling theory accurately describes the dependence of the magnitude of these fluctuations on 
time, on the inter-particle distance and on the size ratio of guest and host particles. The 
consequences in particular for the dynamics of colloidal mixtures are discussed. 



Introduction. - A useful tool in the equilibrium theory of colloidal mixtures is the 'in- 
tegrating out' of the degrees of freedom of one of the components jjj. The resulting renor- 
malisation of the interaction potential between the remaining particles in what has become 
an effective one-particle description, is often referred to as the depletion potential. Depletion 
interactions play a central role in the current understanding of the phase behaviour of disper- 
sions containing colloids of different size and/or shape and have a remarkably elegant physical 
interpretation, although strictly accurate only in the dilute limit. This interpretation hinges 
on the view that an imbalanced osmotic pressure pushes two test particles towards each other 
when they approach within a correlation length set by the other species 0] . 

Depiction potentials are often treated as true, instantaneous potentials even when consid- 
ering time-dependent phenomena such as the kinetics of demixing and gelation |^-^] . This is 
problematical as depletion interactions need time to build up through many collisions between 
a pair of test particles and the particles of the other species. This is (also) the reason why both 
in experiments H, 0] and in computer simulations Q the (mean) depletion potential is only 
obtained after extensive averaging. Here we show by means of scaling theory and Brownian 
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dynamics simulations that the static depletion interaction potential is only approached alge- 
braically in time. Clearly, this must have an impact upon processes such as phase-separation 
kinetics and gelation. 

In this letter we first present the scaling theory of the decay of the fluctuations in the 
depletion force. Next, we outline the simulation technique used and present the results of our 
simulations. We end the paper with concluding remarks. 

Scaling theory. - Arguably, diffusive processes are at the root of the buildup of depletion 
forces, allowing us to estimate the kinetics thereof by fairly simple scaling arguments. For the 
sake of simplicity we consider the situation of two parallel plates in a dispersion of ghost-like 
particles, which do not interact amongst themselves but only with the plates through volume 
exclusion. Our arguments carry over straightforwardly to the more interestingly case of pairs 
of spheres, the results of which we give at the end of this section. 

Considered are two parallel plates of linear dimension a\ separated by a distance R im- 
mersed in a dispersion of (ghost) particles of diameter 02 at a number density p2- The presence 
of the ghost spheres induce a net interaction as they are excluded from the gap between the 
plates when R is smaller than 02 • According to the Asakura-Oosawa depletion theory, JI(J the 
mean interaction force between the plates is given by (F) = Ha 2 for all separations R < 02, 
where IT = k B T P2 is the osmotic pressure of the host dispersion and fc^T the thermal energy. 
For R > cr 2 , (F) = 0. 

The mean number of ghost colloids N 2 (t) contributing to the net force (F) between the 
plates is equal to their mean density times the plate area times a diffusion length, and in- 
creases with time proportionally to p<i<J\\l 'D^t with D2 the diffusivity of species 2 and t the 
time. Assuming the ideal host dispersion to be in local equilibrium, it follows from standard 
statistical mechanics that the mean-square fluctuation of this number must also be of the 
order N%(t). As a consequence, the mean-square fluctuation (AF 2 ) of the depletion force 
should scale as p2k B To\ / \fD$i. Using the Stokes-Einstein relation D2 ~ k B T '/ ' p(J2 for the 
diffusivity of colloid species 2, with p the viscosity of the solvent, we thus obtain 

<AF 2 > ~ p 2 {k B Tf'^lal' 2 p l IH- l l 2 (1) 

for R < 0-2, showing an algebraic decay with time as advertised in the introduction. This is due 
to our assumption of the diffusive nature of the coupling of the plates to the host dispersion. 
On account of the fluctuation-dissipation theorem, eq. (Q) implies that the buildup of the 
depletion forces must also be algebraic in time. From eq. (Q) we furthermore conclude that 
for t t = G2PI 1 p\(k B T') z a\ the depletion interaction between the plates has saturated to its 
mean-field value. The plates act as a microscopic pressure gauge, which equilibrates faster, 
i.e., has a smaller crossover time r, the larger the plates, the higher the temperature, the 
higher the density of the ghost particles, et cetera, as one would in fact expect. 

Our arguments carry over naturally to physically more interesting geometries such as 
pairs of large spheres immersed in dispersions of (ghost) spheres of different size, although 
that the available collision area and the diffusion volume are no longer trivial to evaluate, 
because they now depend on the size ratio of the spheres as well as on their separation (see 
fig. |l|). Here, we simply reproduce without proof for that particular case the final results for 
the mean and the mean-square fluctuation of the depletion force. The former reads (F) = 
7rj 2 (1 + g) 2 [l - R 2 /(rj(l + q) 2 ]p 2 k B T for a x < R < <J X + a 2 and (F) = for R > a t + a 2 , 
where q = (Jija\ is the size ratio. The latter can be shown to obey 

/A/? 2x (f + q)[a 2 (l + q) 2 [l + R 2 /a 2 (l + qf]k B Tf 'p 2 

^ / ~ l(t)[(l + q) + i?/<7i][3<7?(l + q) 2 /A + 3^(1 + q)l{t)/2 + l{tf\ K 1 
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for o\ < R < <j\ +<T2, where l(t) = cyj D^t is a diffusion length. Here, c is an unknown constant 
of proportionality that we fix by fitting to the results of the Brownian dynamics simulations. 
Equation ^| holds also for R > o\ + oi provided we set R = o\ + <72 ; the magnitude of the 
fluctuations in the depiction force level off to that of the net zero force on a single sphere, 
if the separation of the spheres is larger than the diameter of the ghost spheres. Notice the 
non-obvious dependence of (AF 2 ) on time t, distance R and the size ratio q, which we next 
test against the results of Brownian dynamics simulations discussed next. 

Brownian dynamics simulations. - Brownian dynamics simulations were performed using 
the standard scheme developed by Ermak |ll|,^2|. We fixed two test colloids at a distance R in 
a simulation box of volume V, and inserted N2 ghost particles in the same box. The number 
of ghost colloids in the box is determined by the chemical potential in a hypothetical reservoir 
containing only ghost particles. Depending on the size ratio and the density, between 50 and 
f 700 small ghost particles were present in the box. The ghost particles interact with the test 
particles through a hard-core-like steep repulsive potential that scales as e((l + q)ai/r) 50 , 
where e is the interaction strength and r is the distance between the centers of mass of a test 
colloid and a host particle. The force on the particles of species 1 is rendered dimensionless 
by multiplying with er 2 /e. 

Simulation time was divided in a sequence of intervals of width t. For each interval the 
net force f(t) along the vector connecting the centers of mass of the pair of test particles was 
calculated. The variance of the distribution of forces we define as the fluctuation in the force. 
The mean depletion force is simply obtained as (F) = (f(t)), while (AF 2 ) = (f{t) 2 ) - (f(t)) 2 . 
For small values of t the distribution of forces is highly non-Gaussian, very broad and with a 
most probable force that is equal to zero. The distribution only becomes Gaussian for large 
enough t. We fit eq. (|J) to the simulation results by adjusting the parameter c. 

Plotted in fig. || is the value of (AF 2 ) at contact, i.e., for R = o~\, as a function of t for 
q = 0.5, ksT/e = 1 at a ghost particle density P2&2 = 0.0625. The time t is normalised 
on the Brownian time = Snfia^/kBT of the ghost particles, which is the time required 
for a small particle to diffuse over its diameter. Also indicated is the fitted scaling relation 
cq. (||). For small t, the diffusion length is small on the scale of the linear dimension of the 
test particles, and therefore approaches the result for the flat plate discussed earlier. For 
large times the scaling crosses over to i~ 3 / 2 , approaching the time dependence of a radially 
symmetric diffusion volume. Interestingly, best fits were obtained setting C2 = 1 although one 
would expect on basis of the quasi one-dimensional nature of the problem a value of \/2. Note 
that the description is quantitative over 5 decades in time. 

Also in the case of fig. || is the square-mean depletion force (F) 2 . Clearly, (AF 2 ) < (F) 2 
only if t ^> -nf, i.e., for times very much larger than the diffusion time of the host particles. 
This means that for the case q = 0.5, the mean depletion force cannot fully build up at the 
diffusion time of the guest particles. Consequently, the actual depletion force felt by these 
particles must be dominated by fluctuations. 

The dependence of the fluctuations in the depletion force on the size ratio q is given in 
fig. || for a fixed time tDi/a 2 = 0.106 and a fixed number density p 2 o-f = 0.5. Shown are again 
the results for two test colloids at contact. Brownian dynamics simulations and scaling theory 
agree quantitatively and demonstrate that the fluctuations are more prominent the lager the 
size ratio. We were not able to obtain reliable results for q < 0.15 as those simulations required 
an exceedingly small time step to maintain numerical stability. 

Finally we show in fig. [| how the fluctuations vary with the distance between the particles 
for a given size ratio q = 0.5, density P2C1 = 0.5 and tD\ja\ = 0.106. The absolute value 
of the fluctuations in the depletion force increases with the distance between the test colloids 
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and reaches the single particle value for R > o\ + o^- The relative fluctuations (AF 2 ) / (F) 2 
in fact diverge when R — > o~\ + 02 as shown in the inset to fig. |]. This means that fluctuations 
are always important in the tail of the depletion force. Note that for separations close to 
R = (j 1 + (72, the fluctuations are slightly larger than predicted by the scaling theory. This 
discrepancy with scaling theory was not observed for other data sets with different q, t and 
P2, and we have no explanation for it. 

Conclusions. We have shown that the temporal fluctuations in the depletion force 
between two guest particles in a host dispersion of colloids of different size are significant. Our 
results imply that the dynamics of approach of two such colloids must be strongly influenced 
by these fluctuations and the slow buildup of the depletion force. It stands to reason, then, 
that the onset of phase separation or of gelation in mixtures of colloids must be sensitive to 
them. Therefore, one must be wary of treating the depletion potential as an instantaneous 
potential in kinetic studies || ||] . 

The question arises how the fluctuations discussed in this paper can be probed experi- 
mentally One way would be to look at the dynamics of the depletion force using a tweezer 
experiment [fl3| or a TIRM experiment H] , analysing the data in a similar fashion as we have 
done in our simulations. As far as we know this has not been attempted yet. Another way 
would be to look at the escape kinetics of a colloid pushed against a hard wall due to the 
presence of other colloids (similar to the study of which should be strongly influenced 

by a fluctuating depletion force. 

Note added in proof: Upon completion of this work we came across a preprint by Bartolo 
and collaborators jl5fl on fluctuations in medium-induced, Casimir-type forces between parallel 
plates. These Casimir-type forces are due to the coupling of the plates to a fluctuating 
background and are similar but not identical to the depletion forces described in this work, 
because the density fluctuations in our host dispersion are not scale invariant. However, they 
too find that the r.m.s. fluctuation in the effective interaction between them scales as 
as we found in cq. (fil). 

* * * 

This work was funded by EU Grant No. MCFI-1999-00115 [GAV]. We thank Ard Louis, 
Peter Bolhuis, Bob Evans, Paul Bartlett and Gerhard Gompper for useful discussions. 

REFERENCES 

[1] Dijkstra M., van Rou R. and Evans R., Phys Rev. E, 59 (1999) 5744. 
[2] Asakura S. and Oosawa F., J. Pol. Set., 33 (1958) 183. 

[3] Bodnar L, Dhont J.K.G. and Lekkerkerker H.N.W, J. Phys. Chem B, 100 (1996) 19614. 
[4] Verhaegh N.A.M., van Duijneveldt J.S., Dhont J.K.G. and Lekkerkerker H.N.W., 

Physica A, 230 (1996) 409. 
[5] Melrose J.R., Soga K.G. and Ball R.C., J. Chem. Phys., 110 (1999) 2280. 
[6] Poon W.C.K., Renth F., Evans R.M.L., Fairhurst D.J., Cates M.E. and Pusey P.N., 

Phys. Rev. Lett, 83 (1999) 1239. 
[7] Rudhardt D., Bechinger C. and Leiderer P. , Phys. Rev. Lett, 81 (1998) 1330. 
[8] Verma R., Crocker J.C., Lubensky T.C. and Yodh A.G., Phys. Rev. Lett., 81 (1998) 

4004. 

[9] Lowen H., Personal communication. 
[10] Asakura S. and Oosawa F., J. Chem. Phys., 22 (1954) 1255. 



G. A. VLIEGENTHART et al: KINETICS OF DEPLETION INTERACTIONS 5 



[11] We also tested the case of parallel plates against Brownian dynamics simulations which con- 
firmed the scaling relation eq. (|l]). 
[12] Ermak D.L., J. Chem. Phys., 62 (1975) 4189. 

[13] Ohshima Y.N., Sakagami H., Okumoto K., Tokoyada A., Igarashi T., Shintaku K.B., 

Toride S., Sekino H., Kabuto K. and Nishio I., Phys. Rev. Lett., 78 (1997) 3963. 
[14] Kaplan P.D., Faucheux L.P. and Libschaber A. J., Phys. Rev. Lett, 73 (1994) 2793. 
[15] Bartolo D., Ajdari A., Fournier J-B., and Golestanian R., cond-mat/0208095. 



6 



EUROPHYSICS LETTERS 




Fig. 1 - Schematic illustration of the diffusion volume Vd and the collision area A c for two spheres. The 
diffusion volume is defined as a spherical shell of thickness l(t) around a guest particle from which 
the host particles can diffuse towards the surface for collision. The collision area is the available 
area for collisions with a host particle taking into account momentum transfer only along the vector 
connecting the centers of mass of the two guest particles. 
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Fig. 2 - The mean-square fluctuation of the depletion force as a function of the dimensionless time 
t/r2 with T2 the Brownian timescale of the ghost particles. The symbols indicate the simulation 
results and the full curve those of the scaling theory (eq. (^)). The horizontal, dashed line corresponds 
to the average force squared (F) 2 . 
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Fig. 3 - The mean-square fluctuation in the depletion force as a function of the size ratio q for 
pio\ = 0.5, tD\ju\ = 0.106 and R — <T\. Symbols: simulation results, drawn curve: scaling theory 
(eq- (1)). 
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Fig. 4 - The mean-square fluctuation in the depletion force as a function of distance R between the 
two test particles forp20i = 0.5, tDi/af — 1 and q = 0.5. Symbols: simulation results, drawn curve: 
scaling theory (eq. (bh) fitted to the data. 



